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Abstract 

We show that the new precise measurements of Cosmic Microwave Back¬ 
ground (CMB) temperature and polarization anisotropies made by the Planck 
satellite significantly improves previous constraints on the cosmic gravita¬ 
tional waves background (CGWB) at frequencies / > 1CT 15 Hz. On scales 
smaller than the horizon at the time of decoupling, primordial gravitational 
waves contribute to the total radiation content of the Universe. Considering 
adiabatic perturbations, CGWB affects temperature and polarization CMB 
power spectra and matter power spectrum in a manner identical to relativis¬ 
tic particles. Considering the latest Planck results we constrain the CGWB 
energy density to Vt^h 2 < 1.7 x 10 -6 at 95% CL. Combining CMB power 
spectra with lensing, BAO and primordial Deuterium abundance observa¬ 
tions, we obtain Q gw h 2 < 1.2 x 10 -6 at 95% CL, improving previous Planck 
bounds by a factor 3 and the recent direct upper limit from the LIGO and 
VIRGO experiments a factor 2. A combined analysis of future satellite mis¬ 
sions as COrE and EUCLID could improve current bound by more than an 
order of magnitude. 


1. Introduction 

Different processes in the early Universe may have generated a primor¬ 
dial gravitational wave background, such as, among others, quantum per¬ 
turbations during inflation, cosmic strings, pre-big-bang theories, etc (for a 
complete review see [lj and references therein). Detecting this cosmological 
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gravitational wave background (hereafter CGWB) provides a unique way to 
probe the primordial Universe and its evolution. 

The CGWB can be measured at low frequencies constraining a possi¬ 
ble tensor-mode contribution to the large-scale temperature and polarization 
fluctuations in the cosmic microwave background radiation (hereafter CMB). 
The most recent constrain on the tensor-to-scalar ratio is the one published 
by the BICEP2/Keck joint analysis [2], i.e. r < 0.07 at 95% CL [3, 4j, which 
corresponds to O gw /x 2 < 10 -14 in the frequency range 10~ 17 — 10 -16 Liz. 

At higher frequencies, in the range 10 -9 — 10 -8 Hz, pulsars can be used as 
natural gravitational wave detectors, e.g. from the last European Pulsar 
Timing Array (EPTA) data release, Lentati et al. constrain the amplitude of 
GW up to Q gw h 2 < 1.1 x 10~ 9 at 2.8 nHz At smaller scales, interferome¬ 
ters such as LIGO |6j and VIRGO [7j are also looking for gravitational wave 
signals. A recent bound, at ~ 10 2 Hz, is f2 gw /i 2 < 2.6 x 10 -6 from the cross 
correlation between LIGO and VIRGO detectors [8]. Most recently, during 
the review process of this paper, the LIGO-VIRGO collaboration detected a 
Gravitational Waves signal from a Binary Black Hole marger [9]. The cur¬ 
rent best limit on stochastic background coming from Binary Black Holes is 
O gw (/ = 25Hz) = 1.1+H x 1(T 9 at 90% CL [IQ]. 

Moreover, at frequencies greater than ~ 10~ 10 Hz, the stochastic back¬ 
ground can be constrained through big-bang nucleosynthesis (BBN). In fact, 
at these frequencies, primordial gravitational waves contribute to the total 
radiation energy density, increasing the expansion rate of the Universe. In 
this scenario, the CGWB behaves as a free-streaming gas of massless par¬ 
ticles. By measuring primordial abundances of light elements is possible 
to constrain the total number of relativistic degrees of freedom and, conse¬ 
quently, the gravitational waves energy density, for scales smaller than the 
comoving horizon at the end of the BBN m 

Straightforwardly, it is possible to constrain the total radiation density 
through the CMB, reaching even smaller frequencies, ~ 1CV 15 Hz, corre¬ 
sponding to the comoving horizon at the decoupling. In particular, if the 
CGWB energy density perturbations are adiabatic, the extra energy contri¬ 
bution due to gravitational waves is indistinguishable from the one due to 
relativistic neutrinos. Therefore, if we fix the relativistic degrees of freedom 
to its standard value, N e g = 3.046, and parametrize all the extra radiation 
as the effective number of gravitational waves degrees of freedom, N gw , it is 
then easy to translate Ag W into a CGWB energy density, as pointed out in 
(Ej. 
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In this paper we update previous constraints on the CGWB energy den¬ 
sity, as those presented in mm, in light of the latest Planck data release 
and we present the bound reachable combining the future satellite missions 
COrE and Euclid. 

2. Cosmological constraints on GW background 

In this section we discuss the datasets used in the analysis and the ob¬ 
tained results. We make use of the latest CMB Temperature and Polarization 
power spectra of the Planck survey [T^i 16] together with the Planck lensing 
likelihood HU. the Baryonic Acoustic Oscillations observations of [18] and 
the most recent primordial Deuterium abundance observation by Cooke et al. 
[19]. Regarding CMB polarization data at large angular scales, we considered 
both the lowP likelihood based on 70 GHz data, described in [16], and the 
new SimLow likelihood based on 100 and 143 GHz HFI channels as described 
in [20] by including an external prior on the optical depth r = 0.055 ± 0.009 
at 68% CL. 

We explore the cosmological parameters space with the July 2015 ver¬ 
sion of the publicly available cosmomc package m We adopt the following 
parametrization for the ACDM model: the baryon and cold dark matter 
densities oob = and c o c = Q c h 2 , the ratio of the sound horizon to the 
angular diameter distance at decoupling Q MC , the re-ionization optical depth 
r, the scalar spectral index ns, and the overall normalization of the spectrum 
As at k — 0.05 Mpc^ 1 . Furthermore we assume adiabatic initial conditions 
and we impose spatial flatness. We fix the relativistic degrees of freedom, 
parametrized as N e g, to its standard value of 3.046 and, for each MCMC sam¬ 
ple, we compute the primordial Helium abundance assuming standard BBN 
[22], using a recent fitting formula based on results from the PArthENoPE 
BBN code [23]. 

As mentioned above, we parametrize the extra relativistic content due to 
the CGWB as N gw , adding it to the total amount of massless neutrinos. To 
translate to the CGWB energy density, we assume that the contribution 
from a single massless particle to a monochromatic CGWB is 5.6 x ICC 6 , 
therefore (following [TJ 12] and what is done in [241]): 



( 1 ) 
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Figure 1: 95% CL upper limits for the CGWB energy density, tt gw h 2 and for the ef¬ 
fective number of gravitational waves degrees of freedom, N gv , for the different analyzed 
datasets. With Planck-SL we refer to the Planck temperature and polarization Likeli¬ 
hood in combination with the new value of r published in [201 . We report here also the 
value quoted by the LIGO-VIRGO collaboration [5J, the previous cosmological constraint 
obtained combining WMAP7 [22] and SPT[30j, and the Planck 2013 constraint [13] . 


The effective number of gravitational waves degrees of freedom, acting as 
extra radiation, is completely degenerate with the extra relativistic degrees 
of freedom. Therefore, our bounds are model dependent, in fact, any pro¬ 
cess, that modifies the standard radiation content of the Universe, will affect 
the lV gw bounds. For example, the presence of extra relativistic particles at 
recombination (e.g. sterile neutrinos [25]) will tighten the bounds up. On 
the other hand, a reheating phase at MeV temperature can produce a value 
of N eS lower than the standard one [26], relaxing the bounds on lV gw . 

To translate the value of f2 gw /i 2 in to a corresponding value of the tensor-to- 
scalar ratio r is not straightforward since the shape of the tensor spectrum 
needs to be specified. For recent and very comprensive discussions on this 
specific topic see mm- 

We test different combinations of data, starting from Planck 2015 CMB 
data alone (both using lowP and SimLow likelihoods) and then adding BAO 
observations, Lensing data and the primordial Deuterium abundance mea¬ 
surement. We show in figure [l] the 95% CL upper limits for the CGWB 
energy density D gw /r 2 . 
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Datasets 

Cl gw h 2 

N 

1 Vgw 

PlanckTTTEEE - lowP 

< 2.1 x 1(T 6 

< 0.37 

PlanckTTTEEE - lowP + BAO 

< i.9 x nr 6 

< 0.34 

PlanckTTTEEE - lowP + BAO + Lensing 

< 1.6 x 1(T 6 

< 0.29 

PlanckTTTEEE - lowP + BAO + Lensing + Deut. 

< 1.2 x 1(T 6 

< 0.22 

PlanckTTTEEE - SimLow 

< 1.7 x 1(T 6 

< 0.31 

PlanckTTTEEE - SimLow + BAO 

< 1.7 x 1CT 6 

< 0.30 

PlanckTTTEEE - SimLow + BAO + Lensing 

< 1.5 x 1(T 6 

< 0.27 

PlanckTTTEEE - SimLow + BAO + Lensing + Deut. 

< 1.2 x 1(T 6 

< 0.22 

LIGO-VIRGO [8] 

< 2.6 x 1CT 6 


COrE 

< 0.50 x 10~ 6 

< 0.089 

COrE + Euclid 

< 0.076 x 10~ 6 

< 0.013 


Table 1: 95% CL upper limits for the CGWB energy density f l gw h 2 and the effective 
number of gravitational waves degrees of freedom (N gw ) for the considered cosmological 
datasets. We report also bounds from LIGO-VIRGO collaboration, in the frequency range 
41.5 — 169.25 Hz [SJ. With lowP we refer to the likelihood published by Planck in 2015 
m while with SimLow we refer to the new value of optical depth published by Planck 
collaboration [20] . 


Ill Table [T] we report the obtained upper limits on £} gw h 2 and for all 
the data combination that we consider. As expected using the new value 
of r from Planck we reach more stringent constrains. On the other hand it 
is worth to notice that, including the Deuterium abundance measurements, 
this difference is practically negligible. 

The Planck-SimLow alone constraint is 50% better than the upper bounds 
estimated by the LIGO-VIRGO collaboration in 2014 [8]. Combining CMB 
power spectra with the Lensing likelihood and the BAO data we obtain a 
slightly more stringent constraint, Q gw h 2 < 1.5 x ICR 6 ; finally, assuming stan¬ 
dard BBN [[22} [23], we add also the primordial Deuterium abundance mea¬ 
surement reaching fl^h 2 < 1.2 x 10 -6 at 95% CL, improving the precedent 
Planck constraint by a factor of about 3 m, the pre-Planck cosmological 
constraint [13] (based on WMAP [29] and South Pole Telescope [30J results) 
by a factor of 6 and the current interferometer measurements by about 2. 

We also verified the stability of our results with respect to assumptions 
on massive neutrinos. By opening m v as extra parameter, we found, as 
expected, an overall relaxation on bounds. Nevertheless for all the data 
combinations we have considered above, the upper limits increase less than 
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10 %, not changing our conclusions. 

It is also interesting to forecast the future sensitivity on Q gw h 2 achiev¬ 
able with future satellite missions such as COrE |311 32] and Euclid |33i]. To 
this end we simulate mock data for the COrE mission using the 5 channels 
in the range 100-220 GHz, following the approach described in [31], assum¬ 
ing perfect foreground removal and ignoring correlations between multipoles. 
Analyzing this dataset we find that the COrE mission will be able to con¬ 
strain the CGWB energy density to D gw /r < 5.0 x 10 -7 at 95% CL. For the 
Euclid mission we use the fisher matrix approach as described in [351 [36] ; we 
then combine the inverse covariance matrix produced by cosmomc for COrE 
with the Euclid fisher matrix obtaining Q^h 2 < 7.6 x 10 -8 at 95% CL. 

3. Conclusions 

We have used the latest Planck data to constrain a possible cosmological 
gravitational wave background at frequencies greater than 10~ 15 Hz. Our 
tighter constraint is Q gw h 2 < 1.2 x 10 -6 at 95% CL, obtained combining 
CMB with BAO, Lensing and primordial Deuterium observations. This re¬ 
sult improves previous cosmological bounds by a factor 3 and the recent 
LIGO-VIRGO direct measurements by 2. 

We also show that with the next generation cosmological satellite mis¬ 
sions (COrE and Euclid) would be possible to shrink the bounds by more 
than one order of magnitude with respect to current limits. The constraints 
presented here are probably not significant for slow roll inflation that pro¬ 
duces essentially scale invariant spectra. Those models are already strongly 
constrained by current large scale bounds on primordial CMB polarization 
B-modes. However, phase transitions, pseudo scalar inflation or other exotic 
mechanisms that produce a CGWB at higher frequencies (see for example 
ra and [38]) can be constrained by the bounds presented here. 
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